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PSRB1931+24 is the first intermittent radio pulsar discovered to date, charac- 
terized by a 0.8 s pulsation which turns on and off quasi-periodically every ~35days, 
with a duty cycle of ~10%. We present here X-ray and optical observations of 
] PSRB1931+24 performed with the Chandra X-ray Observatory and Isaac Newton 

Telescope, respectively. Simultaneous monitoring from the Jodrell Bank Observatory 
showed that this intermittent pulsar was in the radio-on phase during our observa- 
tions. We do not find any X-ray or optical counterpart to PSRB1931+24 , translating 
Q ■ into an upper limit of ~ 2 x 10 31 ergs -1 on the X-ray luminosity , and of g' > 22.6 on 

I the optical magnitude . If the pulsar is isolated, these limits cannot constrain the dim 

c/) . X-ray and optical emission expected for a pulsar of that age (~1.6Myr). We discuss 

the possibility that the quasi-periodic intermittent behavior of PSRB1931+24 is due 
to the presence of a low mass companion star or gaseous planet, tight with the pulsar 
in an eccentric orbit. In order to constrain the parameters of this putative binary 
' system we re-analysed the pulsar radio timing residuals and we found that (if indeed 

, hosted in a binary system), PSRB1931+24 should have a very low mass companion 

and an orbit of low inclination. 
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1 INTRODUCTION 

A long term radio monitoring study of the ~ 813 ms ra- 
dio pulsar PSRB1931+24 (Stokes et al. 1985; Hobbs et 
al. 2004) , revealed the peculiar intermittent behavior of this 
pulsar (Kramer et al. 2006). PSRB1931+24 is (so far) a 
unique system: it shows an active radio emission phase last- 
ing between 5-10 days (radio-on phase, hereafter), which 
suddenly (in less than 10 s) switches off, and the pulsar re- 
mains undetectable for the following 25-35 days (radio-off 
phase, hereafter). This pattern repeats quasi-periodically, 
and has been monitored for the past 7 years (Kramer et 
al. 2006). Another peculiar property of PSRB1931+24 is 
its spin-down behavior. Remarkably, the pulsar rotation 
slows down considerably faster (by about 50%) when the 
pulsar is in the radio-on phase, with a frequency deriva- 
tive changing from i> on = -(16.30 ± 0.04) x 10~ 15 s" 2 to 
v oS = -(10.80 ± 0.02) x 10" 15 s~ 2 across the two phases. 



From the radio observations, typical pulsar character- 
istics were derived (see Tab. 1 in Kramer et al. 2006): an 
estimate of the dipolar magnetic field (B ~ 2.6 x 10 12 G), 
the characteristic age (t c ~ 1.6 Myr), and the dispersion 
measure (DM = 106.03 ± 0.06 cm" 3 pc). The latter gave 
a rough estimate of the pulsar distance of ~ 4.6 kpc using 
the NE2001 model for the free electron distribution in the 
Galaxy (Cordes & Lazio 2002). Kramer et al. (2006) in- 
terpreted the peculiar spin-down properties of this pulsar 
in terms of a transient plasma flow in the magnetosphere, 
causing the quenching and re- ignition of the radio emission. 
In particular, an increase of the magnetospheric current flow 
during the radio-on phase, can provide an additional brak- 
ing torque on the neutron star. From the observed variations 
in the pulsar spin-down, these authors showed that the pul- 
sar magnetospheric electron density was consistent with the 
Goldreich- Julian density (Goldreich & Julian 1969). This 
interpretation can successfully explain the variable neutron 
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star torque, but still fails to explain what causes the changes 
in the magnetospheric plasma flow, especially if in a quasi- 
periodic fashion. To date, this periodicity of the radio-on 
and radio-off recurrence is difficult to explain in any sce- 
nario considering an isolated pulsar. As Kramer et al. (2006) 
pointed out, the short shut-off time of less than 10 seconds, 
and the rather stable pulse profile, rule out possible scenarios 
like precession of the neutron star. On the other hand, the 
long radio-off phase is in contrast with the typical nulling 
timescales of radio pulsars, exceeding it by almost five orders 
of magnitude. 

Cordes & Shannon (2006) studied the possibility that 
PSRB1931+24 might be in a binary system with a large 
asteroid, surrounded by a disk of small asteroids, moving 
in a ~ 40 day eccentric orbit around the pulsar. This sce- 
nario suggests that the interaction between the pulsar mag- 
netosphere and the asteroids can be responsible both for 
the torque change and the intermittent radio activity of 
PSRB1931+24. Recently, Zhang, Gil & Dyks (2007) have 
proposed that intermittent pulsars are old isolated neutron 
stars which entered (or are about to enter) the so called 
"death valley" (Chen & Ruderman 1993), where the po- 
lar cap voltage of the pulsar is not sufficient to power 
the pair production process. However, this scenario still 
fails to account for the quasi-periodic pulsed emission of 
PSRB1931+24. 

Here we report on X-ray (Sj2j and optical (JU obser- 
vations of PSRB1931+24 taken with the Chandra X-ray 
Observatory and the Isaac Newton Telescope. Furthermore, 
in $4] we derive constraints on possible orbital parameters 
and companion mass by studying the radio timing residu- 
als. In Sj5] and we discuss our results and investigate in 
detail the possibility that PSRB1931+24 resides in a bi- 
nary system with a ~35 day orbital period and a low mass 
companion star or a planet. 
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Figure 1. Finding-chart of the i' field of PSRB1931+24. The 
circle is the uncertainty in the pulsar postion as derived from 
radio timing (see Tabl). 

of 7x 10 -15 ergs -1 cm -2 or 1.2 x 10 -14 ergs -1 cm -2 , assum- 
ing a black body (kT = 0.3 keV) or a power law (P = 2.5) 
spectral decomposition, respectively. At 4.6 kpc these fluxes 
translate into a 99% upper limit on the 0.3-10 keV X-ray 
luminosity of 1.7 x 10 31 ergs -1 and 3 x 10 31 ergs -1 , depend- 
ing on the assumed X-ray spectrum. However, note that the 
distance of 4.6 kpc, as inferred from the pulsar DM (Kramer 
et al. 2006) might well have an uncertainty of a factor of 2 
due to a ~30% uncertainty in the Galactic electron density 
model (Cordes & Lazio 2002). Hence, more conservative up- 
per limits on the luminosity would be 4 times larger than 
reported above. 



2 X-RAY OBSERVATION 

The Chandra Advanced CCD Imaging Spectrometer (ACIS) 
observed PSRB1931+24 on 2006 July 20th, for an on-source 
exposure time of ~9.7ks. The target position was placed in 
the standard back-illuminated ACIS S3 aimpoint, using the 
FAINT mode. We corrected the astrometry for any process- 
ing offset and we cleaned the image for hot pixels. Run- 
ning CIAO celldetect and wavedetect toolflj, no X-ray 
sources were detected in the whole ACIS-S3 CCD, while 
four (unrelated) sources were detected in the other CCDs 
at > 5cr confidence level over the background. In partic- 
ular, no photons were detected in a 1" radius around the 
PSRB1931+24 radio position (see Tab. 1 for the RA and 
DEC). We then obtained a 99% upper limit on the source 
count rate of 4.74 x 10 -4 counts -1 (Gehrels 1986). 

Assuming a conservative absorption value of Nh = 
8.3 x 10 21 cm -2 (derived for the pulsar position from Dickey 
& Lockman 1990, hence considering the whole Galactic Nh 
in that direction), and using the PIMMS calculatoiQ, we de- 
rived a 99% upper limit on the 0.3-10 keV unabsorbed flux 

1 for details refer t o |http: / / asc.harvard.edu/ciao J\ 

2 for details refer to http: / /heasarc. gsfc.nasa.gov/Tools/w3pimms. 



3 OPTICAL OBSERVATIONS 

We observed the field of PSRB1931+24 with the Isaac 
Newton Telescope located at the Roque de Los Muchachos 
Observatory in La Palma, for an exposure time of 10 min- 
utes in three optical filters: g', r' and i' (see Fig. [1] for the 
finding-chart). No optical counterpart was detected within 
the radio pulsar position. We derived 5a upper limits on the 
optical emission of PSRB1931+24 of g' > 22.6, r' > 22.4 
and i' > 22.2 magnitudes. We inferred these optical upper 
limits from the magnitudes of the faintest object detected 
at 5a confidence level in the same CCD as the pulsar posi- 
tion. For the astrometry we used an 11' x 11' subsection of 
the g' image where we found 85 stars from the USNO CCD 
Astrograph Calalogue (UCAC2; Zacharias et al. 2004). We 
obtained an astrometric solution, fitting the zero-point po- 
sition, scale and position angle; the final rms residuals were 
0".13 in both in RA and DEC. 

We inferred the reddening in the direction of 
PSRB1931+24 from the Nh, which gave an A v = 4.64 
(Predehl & Schmitt 1995). We then converted this value into 
an estimate of the reddening in the filter we actually used 
(Rieke & Lebofsky 1985; Schlegel et al. 1998): A g , = 4.94, 
html|A r / = 3.68 and A t i = 2.83. Considering a distance of 4.6 kpc, 
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Parameter 


Value 


RA (J2000) 


19:33:37.88(5) 


DEC (J2000) 


24:36:40(1.5) 


Epoch 


52281.296676 


v (Hz) 


1.22896706877(3) 


;> (icr 15 s" 2 ) 


-12.1501(4) 


v (10 -25 s" 3 ) 


-2.0(2) 


DM (cm -3 pc) 


106.03(6) 



Table 1. Parameters derived from the radio timing analysis de- 
scribed in ^4] Errors are at la confidence level. 



we inferred upper limits on the source absolute magnitude 
of: g' > 4.35, r' > 5.2 and i' > 6.1. From our upper lim- 
its of the absolute optical magnitudes, we could derive a 
rough range of stellar types for a putative counterpart. From 
r' > 5.2, which was the most constraining limit, we could 
exclude all giant and super-giant stars, and deriving a spec- 
tral type later than a G8 star. Note that having we assumed 
the absorption value of the whole Galactic Nh'- any (more 
realistic) lower Nh values would result in stellar types much 
later than a G8. 



light-seconds, which is a function of a p — hr M c / (M p +M c ), 
the semi-major axis of the orbii0- 

On the other hand, the relative orbit art, the orbital 
period and masses are related according to Kepler's 3rd law: 

^(^) 3 = T (M P + M C ), 

^orb V C / 

with the masses measured in solar units and T Q = 
GM & /c a = 4.925490947^s. Assuming P olb ~ 35 days and 
M p ~ IAMq we could derive an estimate for the orbital 
inclination angle as a function of the companion mass in 
the following way: for a given companion mass M c we per- 
formed Monte- Carlo simulations drawing possible values for 
E, uo and e from uniform distributions over [0, 2n] and [0, 1], 
respectively. Note that the distribution for the eccentricity 
is not expected to be uniform (i.e. it is likely to be skewed to 
small eccentricities) but for our purposes this assumption is 
sufficient as we were mostly interested in upper bounds for 
the inclination angle. For each companion mass, one million 
Monte-Carlo runs were performed and the median and 95% 
confidence limits were computed (see Fig. [3] the grey shad- 
owed region is the allowed region at 95% confidence level). 
Hence, for reasonable companion masses the tight limit on a 
detectable periodicity in the timing residuals implies rather 
small orbital inclination angles (see also fjGj) . 



4 LIMITS ON THE ORBITAL PARAMETERS 
FROM THE RADIO TIMING 

If PSRB1931+24 happens to be in a binary system (e.g. 
as investigated in detail in §[S]), the motion of the pulsar 
around the system's centre of mass should leave a periodic 
signature in the remaining timing residuals. We searched for 
this signature in radio timing data spanning almost 7 years 
(Kramer et al. 2006). We fit the radio timing residuals with 
one spin frequency derivative and a second derivative (values 
are reported in Tab. 1). The latter was used to remove a cu- 
bic term due to long term timing noise. In order to improve 
the position information, we removed the remaining timing 
noise by using the fit-wave method described by Hobbs et 
al. (2004), with a minimum period of 2 years. The pre- fit- 
wave residuals were then analysed for periodicities caused 
by low-mass companions (with masses smaller then the pul- 
sar mass, M c « M p ), as described by Freire et al. (2003). 
As Fig. 2 shows, we detected, unsurprisingly, a faint signal 
at a period of ~ 37 days. However, we believe this signal is 
consistent with a periodicity caused by the sampling of time- 
of-arrivals forced upon us by the quasi-periodic on-phases of 
the pulsar. We do not detect any other significant signal. 

For a more general companion mass we also performed 
the following analysis of the timing residuals. The peak- 
to-peak amplitude of the remaining timing residuals was 
Ai res ~2 ms. From this we placed limits on the correspond- 
ing orbital parameters by interpreting At TCS as caused by a 
"Roemer delay", i.e. the light-travel time across the orbit 
(e.g. Lorimer & Kramer 2005). We then set: 



= x (cos E — e) sin uo + sin E\f\ 



equal to Aires. Here, E is the eccentric anomaly, e the eccen- 
tricity and lo the longitude of the periastron. Furthermore, 
x = a p sini/c is the projected semi-major axis measured in 



5 INVESTIGATING THE POSSIBLE BINARY 
NATURE OF PSRB1931+24 

In this section we investigate the possibility that the pe- 
culiar characteristics of PSRB1931+24 can be explained 
by a low mass star or a gaseous planet moving in an ec- 
centric orbit around the pulsar (see also §13]). This causes a 
quasi-periodic intermittent behavior, as well as the change in 
the spin-down characteristics of the pulsar, in two different 
ways: 1) providing at the periastron passage the additional 
plasma needed in the Kramer et al. (2006) model, or 2) ac- 
creting material on the pulsar magnetosphere and causing 
the transition between the radio pulsar and the propeller 
regime (see for details Illarionov & Sunyaev 1975; Stella et 
al. 1986, 1994). In this section we discuss only this second 
possibility. The physics involved in the first possibility has 
been thoroughly investigated by Kramer et al. (2006), al- 
though without considering the wind of a companion as the 
source of the additional magnetospheric plasma during the 
radio-on phase (see also §[S]). 

We define here three important PSRB1931+24 radii 
which will be used in the following. 
The magnetospheric radius: 



1.78 x WM[ 



■2/7 



the corotation radius: 



0.12 x 10 8 M^P^ 3 



2.25 x 10 s cm 



and the light cylinder radius: 



R lc = 0.46 x 10 8 Pio ^ 3.74 x 10 9 cm . 

3 i is the orbital inclination angle, c is the speed of light, and 
is the size of the relative orbit. 
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Figure 2. Lomb-Scale spectrum of PSRB1931+24. Dashed 
curve indicates a 99.99% confidence level. A periodicity is de- 
tected at ~ 37 days, but see text for details. 

These three radii represent, respectively: i) the balance of 
the gravitational force and the ram pressure of the infalling 
material, ii) the centrifugal barrier for the infalling material 
due to the pulsar rotation, and iii) the place where the pulsar 
magnetic field lines open such that their tangential velocity 
do not exceed the speed of light (Illarionov & Sunyaev 1975; 
Ruderman & Sutherland 1975). 

We define: Bg — B ns /10 9 G is the neutron star mag- 
netic field, Pio is the spin period in units of 10 ms, M15 = 
M/10 15 g s _1 is the mass inflow rate toward the neutron 
star, and M 1A = M ns /1.4M Q and Re = R ns /W 6 cm are the 
mass and the radius of the neutron star, hereafter assumed 
Mi. 4 = Re = 1. Note that during the orbital motion, the 
only variable radius is i? m oc M^ 5 . 



5.1 Radio pulsar is on 

If the mass inflow rate from the putative companion star is 
very low, as might happen at apastron of an eccentric orbit, 
the magnetosphere of the neutron star might be larger than 
the light cylinder and corotation radii (i? m > R\ c > R COI ). 

As a result, the centrifugal barrier is closed, which 
means that incoming material from the companion star is 
prevented from falling towards the neutron star magneto- 
sphere or surface. In this case the neutron star might behave 
as a radio pulsar. The pulsar radiation pressure dominates 
over the ram pressure of the inflowing material, preventing 
the matter to penetrate towards the neutron star (Illarionov 
& Sunyaev 1975; Davis & Pringle 1981; Stella et al. 1994; 
Campana et al. 1998). 

This occurs when the mass inflow rate (M) is smaller 
than the limiting value: 

M 15 ,radio-off = 5.4 x lfT 2 BlP{ 7/2 M-l /2 Rl ~ 0.075 , 

above which the radio pulsations would quench, even as 
rapidly as ~10s. 

It is worth noting that when M ~ Mi5 jra dio-off , the 



radio pulsar emission is not in equilibrium (Illarionov & 
Sunyaev 1975; Shaham & Tavani 1991), depending on the 
stability of the mass inflow rate; sporadic variability of the 
radio-on duration of PSRB1931+24 might take place, as 
indeed observed (Kramer et al. 2006). Furthermore, during 
the radio-on phase, variations in the pulsar DM as a func- 
tion of the orbital phase are expected. Unfortunately the 
observations were performed at one frequency, making im- 
possible to put a meaningful limit on the DM variability (see 
also Kramer et al. 2006). 



5.2 Radio pulsar is off 

If the mass inflow rate starts to increase, e.g. approaching 
periastron, we expect a correlated decrease of the magneto- 
spheric radius, which eventually becomes smaller than the 
light cylinder radius (Ri c > i? m > R CO i)- This happens when 
the mass inflow rate towards the neutron star becomes larger 
than the limiting value Mi 5iI . a dio-off ; at this point the pulsar 
radiation pressure is overcome by the ram pressure of the in- 
falling material quenching the radio pulsar mechanism. Dif- 
ferent regimes are then allowed at this stage. If i? m remains 
larger than R COI , which means Mi 5 , ra dio-off < M < Mis,acc, 
with 

Mi 5 ,acc = 5.97 Blp- 7/3 M-2 /:i R b 6 ~ 1.4 x 10 3 , 

then the magnetosphere of the neutron star still rotates in 
a super-Keplerian motion, and the inflowing material might 
either accumulate outside the magnetospheric boundary or 
be swept away by the magnetospheric drag: this is called 
the "propeller" regime (Pringle & Rees, 1972; Illaroniov & 
Suniaev 1975; Davies & Pringle 1981; Wang & Robertson 
1983; Stella, White & Rosner 1986). 

Given the upper limits we derived for the possible com- 
panion star (see §[3]and §|4|, there is only a very small chance 
that the M, due to the wind loss of the companion, over- 
comes Mi5 iacc . We will then consider hereafter only the pos- 
sibility that the radio-off phase is driven by the propeller 
regime, excluding the surface accretion scenario. 

When the radio pulsations are quenched, the spin down 
behavior of the pulsar is not driven anymore by the mag- 
netic dipolar loss as it was before. What happens to the 
pulsar spin-down during the propeller regime is still rather 
controversial, and requires detailed hydrodynamic simula- 
tions (Romanova et al. 2003). Depending on which kind of 
instability and shock takes place on the pulsar magneto- 
sphere, the spin-down rate might either increase further or 
be reduced by the angular momentum transferred by the 
infalling material to the magnetosphere. 

In the PSR B1931+24 case, it is clear that the infalling 
material should provide a certain rotational energy in or- 
der to make the pulsar slow down less efficiently during the 
radio-off phase, changing the spin-down of the neutron star 
by Av = i> on - i> oS = -5.5 ± 0.4 x 10" 15 Hzs -1 , which 
converted in energy corresponds to AE ~ 4ir 2 IuAv ~ 
2.6 x 10 30 ergs _1 . 

When the infalling matter reaches the neutron star mag- 
netosphere, the source is expected to emit in the X-ray band 
with a minimum luminosity of: 

Lradio-ofi ■ = 7 x lQ zl Blp- s,2 Ml'2Rl m 0.12 x 10 31 ergs~\ 
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Figure 3. Limits on the orbital inclination angle as a function 
of the companion mass as derived from the pulsar radio timing. 
For each companion mass, one million Monte-Carlo runs were 
performed, the median and 95% confidence limits were computed. 
The grey shadowed region is the allowed region at 95% confidence 
level. 



and a maximum luminosity related to the maximum wind 
loss we expect from the inferred companion type. For a 
stellar type later than a G8, we expect the X-ray lumi- 
nosity of the system in the radio-off phase to be between 
1.2 x 10 30 - 1.6 x 10 32 ergs~\ 

After the periastron passage, the low-mass companion 
star moves towards apastron, the magnetospheric boundary 
expands because of the decreasing mass inflow rate, the cen- 
trifugal barrier closes and the radio pulsar mechanism can 
resume again. 



6 DISCUSSION 

The upper limits we derived from these X-ray and optical 
observations of PSRB1931+24 during the radio-on phase, 
are not deep enough to detect the typical X-ray luminosity 
(< 10 29 erg s _1 ; Becker & Triimper 1997) and optical mag- 
nitude (V>28; Mignani et al. 2007) of an isolated pulsar of 
~1.6Myr, as PSRB1931+24. 

On the other hand, if the pulsar is hosted in a binary 
system, we derived tight constraints on the putative com- 
panion and its orbital parameters. We constrained the com- 
panion to be a low mass star later than a G8 type or a 
gaseous planet. For an orbital inclination of 0.1 degree the 
companion would have a mass M c <0.01 Mq, while with 
a smaller inclination of 0.01 degree, the companion mass is 
allowed to be M c <0.1 Mq (all at 95% confidence level; see 
Fig. 3). 

In the binary scenario (see §[5| , the pulsar intermittence 
can be explained by the transition between the radio pulsar 
regime and the propeller regime during the orbital motion. 
The X-ray upper limits we derived are unfortunately rel- 



ative to the radio pulsar dominated regime (note that the 
source was radio-on during the observations), and not to 
the propeller regime when the source can be bright in the 
X-ray due to accretion of matter onto the pulsar magne- 
tosphere. Our Chandra observation was in fact aimed at 
detecting a possible X-ray emission from PSRB1931+24 
during its radio-off phase. Unfortunately, the relatively un- 
certain duration of both phases (the radio-on phase rarely 
might happen to be longer than the average duration), and 
the relatively large error (~ a few days) on the putative pe- 
riastron passage, made the observation be performed during 
the radio-on phase of the pulsar. 

Note that the insufficient flux detection limit of cur- 
rent and past X-ray surveys and monitoring programs (e.g. 
ROSAT and the RXTE/ASM), do not make the non detec- 
tion of an X-ray counterpart in previous X-ray surveys at 
all constraining. 

Although very successful in explaining the intermit- 
tence, the quasi-periodicity and the torque variations of 
PSRB1931+24 , this binary scenario still presents several 
unsolved issues. One issue of this scenario, and to some ex- 
tent also to the "asteroid" model (Cordes & Shannon 2006), 
is the very small orbital inclinations required by the radio 
timing studies (§[4] and Fig. 3). For a random distribution 
of inclination angles, the probability of observing a system 
at an angle less than a value io, is p(i < in) = 1 — cosio- 
Considering the favorable case of an angle of 15° (a gaseous 
planet; see Fig. 3), the probability of observing such a sys- 
tem is p(i < 15°) = 3.4%. A chance of 3.4% is not very high, 
although still worth to be considered. 

Another issue concerns the accretion rate from the com- 
panion. Can the accretion from the low wind loss of such a 
small mass star or gaseous planet be high enough to switch 
the pulsar from the radio emitting regime to the propeller 
regime? Typical wind rates for e.g. a K0 star are insuffi- 
cient to produce the limiting value of Mi 5ira( jio-off by ~2 
orders of magnitude. This problem might be partially (but 
not totally) alleviated by taking into account the irradiation 
process (Podsiadlowski 1991; DAntona & Ergma 1993). In 
particular, a low mass star irradiated by the radio pulsar is 
expected to expand even if the radiation bath is not partic- 
ularly extreme, and the wind of the star to increase substan- 
tially. However, the irradiation process cannot be dominant 
on the stellar wind from a G8 to K0 star, mainly because of 
the relatively low rotational energy of PSRB1931+24 with 
respect to the surface temperature of the companion star 
(Li rr — fE Iot (Rs/2a p ) 2 , where R s is the companion star 
radius). For lower mass stars or considering a gaseous planet 
orbiting around the pulsar, this process might have instead 
a substantial role, and might produce strong winds towards 
the pulsar while at periastron. 

If the weak stellar wind from the low-mass companion 
star is insufficient to drive the radio pulsar in the propeller 
regime, it is instead strong enough to provide the amount 
of plasma needed to produce the spin-down change in the 
scenario proposed by Kramer et al. (2006) . Note that in this 
picture, at the periastron passage we expect the radio-on 
phase, while in the propeller scenario it is expected at the 
apastron of the orbit. However, while the propeller scenario 
naturally explains the quenching of the radio emission, in 
this picture the radio quenching still remains puzzling: we 
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would in fact expect to see the pulsar at all the time, al- 
though with two different spin-down behavior. 
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